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Introduction: At present, visual inspection is the only method for comparing the large-scale 
morphologies in the distribution of galaxies to those in model universes generated by N-body simulations. 
To remedy the situation, we have developed a set of three structure functions (Si, S 2 , S 3 ) that quantify 
the degree of large-scale prolateness, oblateness and sphericity /uniformity of a 3-D particle distribution 
(c./. Babul & Starkman 1992, preprint) and have applied them to a volume-limited (< 4000 kms -1 ) sample 
of 699 IRAS galaxies with f^o > 1.2 Jy (Fisher 1992, Ph.D. Thesis, U.C. Berkeley). 

To determine the structure functions, we randomly select 500 galaxies as origins of spherical windows 
of radius Rw , locate the centroid of the galaxies in the window (assuming all galaxies have equal mass) and 
then, compute the principal moments of inertia (/ 1 , I 2 , I 3 ) about the centroid. Each Si is a function of I 2 /I 1 
and J 3//1 . Si, S 2 and S 3 tend to unity for highly prolate, oblate and uniform distributions, respectively and 
tend to zero otherwise. The resulting 500 values of S,- at each scale Rw are used to construct a histogram. 

IRAS Sample: In analyzing the IRAS sample, we restrict ourselves to windows with radius 
1100 kms -1 < Rw < 1700 kms -1 . Smaller windows contain too few galaxies while larger windows are 
rejected due to the scarcity of window origins for which the window lies entirely within the sample volume. 
A set of histograms for Si , S 2 and S 3 in windows of radius 1300 kms -1 is shown in panel 1 of Figure 1. The 
structures in the IRAS galaxy distribution do not appear to be highly prolate or oblate. 

Model Distributions: We compare the IRAS results to analyses of eight galaxy distributions 
discussed by Weinberg & Cole (1991, preprint). The eight distributions correspond to particles identified 
as galaxies according to biased and unbiased schemes from simulations of an fl = 1 universe characterized 
by the initial 1-point probability distribution, P( 6 ), which is either Gaussian, or has a long tail of positive 
fluctuations (skew-positive), a long tail of negative density contrasts (skew-negative), or both negative and 
positive tails (broad). The non-Gaussian initial conditions were generated by a local transformation of the 
Gaussian field with power spectrum P(k) oc k~ l . In order to ensure that 5, for the models are comparable 
to those for the IRAS galaxies, we analyze four subsets (spheres of radius 4000 kms -1 ) of each of the 
simulation volumes, each containing 400 — 1000 galaxies. Placing the observer at the centre of the spheres, 
we transform the simulations from real space to redshift space. For each of the four volumes, we compute 
the Si histograms at Rw = 1100, 1300, 1500, and 1700 kms -1 . For each 5,-, the four histograms at a given 
Rw are then averaged to obtain the mean “model” histogram. Panels 2 and 3 of Figure 1 show the mean 
histograms for the unbiased Gaussian and the unbiased skew-positive models (Rw = 1300 kms -1 ). 

Comparison: As a prelude to comparing the IRAS histograms to those for the models, we interpret 
the mean model histograms (for each 5,- and at each Rw) as the underlying probability distributions and 
compute the likelihood (using the multinomial likelihood function) of obtaining the corresponding IRAS 
histogram. We also compute the likelihood for obtaining the histograms for each of the four spherical sub- 
volumes from the mean. The results for the two least unlikely models (unbiased skew-positive and unbiased 
Gaussian) are presented in Figure 2. The skew-positive model is preferred; its likelihood is the largest of 
any of the models and it is the only one for which the IRAS results are comparable to the likelihoods of 
individual volumes for all three Si and at all Rw- The results of our likelihood analysis are evident in the 
histograms themselves. For example, in the Gaussian model, the histogram for Si (Rw = 1300 kms -1 ) 
has a tail extending to Si = 0.3 while both the skew-positive and IRAS histograms (compare panels in 
Figure 1) are narrower. The differences between the 53 histograms are even more pronounced, with the 
Gaussian histogram being noticeably broader than either the IRAS or the skew-positive. The histograms 
for the latter two are similar. As an aside, we note that four out of five colleagues, when asked to visually 
compare unlabelled particle plots of the model galaxies to those for IRAS galaxies, also selected the unbiased 
skew-positive model as the most likely match. 

Conclusions A preliminary comparison of the large-scale morphology of IRAS galaxies with those 
present in the eight simulated galaxy distribution suggests that the distribution of IRAS galaxies is most like 
the distribution of galaxies in an unbiased, evolved skew-positive model. This result suggests that either the 
initial conditions that gave rise to the IRAS galaxies are non-Gaussian or the nature of biasing associated 
with their formation gives such an impression. 
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Fig. 1. Histograms for the normalized frequency distribution of 5,- for IRAS galaxies (left column of 
panels) as well as for the galaxy distribution (unbiased) from simulations with Gaussian (middle column) 
and ppsitively skewed (right column) initial 1-point probability distribution. P(6). 



Fig. 2. The value of the likelihood function for obtaining the IRAS S’j histograms at IZu, = 1100, 1300, 1500, 
and 1700 krns" 1 from the histograms for the skew-positive model (solid line) and for the Gaussian model 
(dashed line). The points show to the likelihood of obtaining the histograms for the four sub-volumes for 
each model from the mean histogram. The open circles correspond to the skew-positive model and should 
be compared to the solid line; the crosses correspond to the Gaussian simulation and should be compared 
to the dashed line. 
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